Y 
© 


Y 
© 


wo g 
a ô 


Pulse-shaping method for real-time neutron/gamma discrimination at low sampling rates“ 
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The Nal:TI scintillator is an innovative material for dual-gamma-ray and neutron detection with a low °Li 
concentration. To achieve real-time n/y discrimination, a zero-crossing time comparison algorithm based on 
trapezoidal pulse shaping was developed. The algorithm can operate efficiently at low sampling rates and was 
implemented on a single-probe portable digital n/y discriminator based on a feld-programmable gate array 
(FPGA). The discriminator and Nal:TI,°Li detector were tested in a neutron-gamma mixed field produced by 
an 741 Am-Be neutron source to evaluate the performance of the algorithm. The Figure of Merits (FoM) was 
measured as 2.88 at a sampling rate of 50 MHz, indicating that the discriminator with its embedded algorithm 
has a promising n/y discrimination capability. Efficient discrimination at sampling rates of 40 and 25 MHz 
demonstrates that the capability of this method is not limited by low sampling rates. 
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I. INTRODUCTION 


Efficient neutron detection is in high demand in numerous 
fields, such as radiation monitoring[1], neutron imaging[2], 
and national security[3]. *He gaseous detectors are widely 
used in neutron detection systems. Although *He position- 
sensitive detectors show excellent neutron-gamma discrimi- 
nation ability, reasonable spatial resolution (1-10 mm), and 
high count-rate capability[4], the high cost and shortage of 
3He limit their application. The world supply of ?He used 
to come entirely from the decay of tritium originating from 
nuclear weapons programs in the U.S. and Russia[5]. The 
mainstream choice of neutron detectors has shifted to scintil- 
lator detectors, including organic and inorganic scintillators. 
Because neutron emission is often accompanied by gamma 
rays, n/y discrimination is a critical technique for neutron de- 
tection. 

Scintillator detectors not only provide a high light yield 
and excellent energy resolution but also emit discrimina- 
tive neutron pulse signals that can be extracted by pulse 
shape discrimination (PSD) methods[6]. Most PSD meth- 
ods have been tested with plastic scintillators, such as EJ- 
276 and EJ-309, and most tested FoMs are better than 1.0[7- 
10], which means that they can extract most neutron signals 
from mixed neutron-gamma radiation field. Plastic scintil- 
lators are preferred over liquid scintillators because of their 
attractive features, including low cost, self-containment, and 
ease of processing [10]. Among the anthracene, stilbene, and 
p-terpheny] scintillators, anthracene scintillators exhibit the 
highest light yield and PSD performance[11]. Inorganic scin- 
tillators have been developed using °Li(n, «)?H reaction. The 
reasonable discrimination performance of elpasolite scintilla- 
tors, such as CLYC, has been verified.[12, 13]. As the latest 
inorganic scintillator, Nal:T1,°Li is capable of efficient dual 
detection of thermal neutrons and gamma rays by doping with 
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SLi while maintaining the original NaI:T1 scintillation prop- 
erties. Moreover, Nal:T1,°Li scintillators can be significantly 
larger than CLYCs[6]. 


In recent years, many time- and frequency-domain meth- 
ods have been proposed and compared based on different fea- 
tures of the output signals from scintillators[7]. For example, 
the charge-comparison method integrates pulse signals and 
calculates the charge ratio in the time domain[7], whereas 
the discrete wavelet transform method calculates a scaling 
function that represents the energy density in the frequency 
domain[14]. Additionally, intelligent techniques such as neu- 
ral networks have been applied to PSD [9, 15]. However, 
many matrix operations may be required in the frequency do- 
main and artificial intelligence methods, which are limited 
to offline computing[15]. Such tedious and complex algo- 
rithms require long identification times; thus, they are not su- 
perior to fast real-time analysis. A confusing problem with 
time-domain methods is that their PSD performance is not as 
satisfactory as that of intelligence methods. One approach 
to solving this problem involves the use of pulse shaping, 
which can be applied to alter the features of the output sig- 
nals from scintillators. A larger difference in signal shapes 
leads to better PSD performance. The pulse-shaping tech- 
nique can decompose and synthesize signals in real time in an 
FPGA[16, 17] and potentially simplify PSD methods while 
maintaining good PSD performance. 


Several digitizer boards are used to sample signals in most 
traditional PSD methods; however, a host computer is re- 
quired for data processing[18—20]. A possible constraint 
for developing a real-time discriminator without a computer 
is the adaptability of the PSD algorithms. Additionally, to 
achieve better PSD performance, some algorithms have high 
hardware requirements [19-24], such as high sampling rates 
and large FPGA logic resources. However, the device size, 
power consumption, and cost have also been considered in 
some applications[25]. This paper proposes a time-domain 
algorithm based on a trapezoidal shaper to achieve real-time 
n/y discrimination and demonstrates the feasibility of deploy- 
ing this algorithm in a portable FPGA-based discriminator. 
The PSD performance was evaluated using a Nal:TI°Li scin- 
tillator at low sampling rates. 
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Il. METHOD 


A. Hardware design 


A block schematic of the n/y discriminator is shown in 
Fig. l(a). The n/y discriminator receives signals from the 
Nal:TI,°Li detector and uploads the latest spectra to a com- 
puter every two seconds. Most signal processing is performed 
in FPGA so the analog circuit only needs to satisfy the sig- 
nal amplification and noise reduction requirements. After 
the neutron and gamma pulse signals are produced by the 
preamplifier, they must be further amplified and adjusted in 
the analog circuits of the discriminator. AD829 achieves its 
230 V/s uncompensated slew rate and 750 MHz gain band- 
width while requiring only 5 mA of current from power sup- 
plies, which is suitable for building a main amplifier. With 
a low input voltage noise of 1.7 nV/VHz, AD829 serves 
as an input buffer for an analog-to-digital converter (ADC). 
The amplified analog signals are then digitalized in a high- 
speed ADC (AD9236, Analog Device), which uses a multi- 
stage differential pipelined architecture to provide 12-bit res- 
olution at up to 80 Msps. AD9236 operates from a single 
3.3 V power supply and has a low 366 mW power consump- 
tion at 80 Msps. The sample-to-hold amplifier (SHA) input 
maintains excellent performance for input frequencies up to 
100 MHz and is configured as single-ended in this appli- 
cation. A single-ended clock input was used to control all 
the internal conversion cycles provided by FPGA. A clock 
duty cycle stabilizer (DCS) was used to re-time the non- 
sampling edge and provide an internal clock signal with a 
nominal 50% duty cycle. This allows for a wide range of 
clock input duty cycles without affecting the performance of 
AD9236. Hence, additional single-ended-to-differential cir- 
cuits or clock sources are not required. Using an internal 
reference voltage, the circuit design can be simplified. Fi- 
nally, the digitalized signals were delivered to the Spartan-7 
FPGA (XC7S25-1CSGA324C, Xilinx), where a critical dis- 
crimination algorithm was implemented. The FPGA must 
perform three main tasks: provide the clock source for the 
ADC, receive and process the digitized signal from the ADC, 
and upload the energy spectra. Digital circuits require volt- 
ages of 1.0V, 1.8V, and 3.3V. These low voltages were gener- 
ated using low-dropout (LDO) linear regulators (TLV74110, 
TLV74118, TLV75533, TI). These LDO regulators have fixed 
output voltages and their peripheral circuits require only ca- 
pacitors for filtering. 

The PSD algorithm was implemented using several mod- 
ules on FPGA. The digital antinoise filter and signal prepro- 
cessing module use FIFO to cache the data from the ADC. 
Furthermore, they filter the high-frequency noise and adjust 
the pulse shape when severe pulse undershoot occurs. In- 
stead of being sent directly to the discrimination module, the 
filtered signals are first processed using a trapezoidal pulse 
shaper. The trapezoidal pulse shaper is conducive to mag- 
nifying the difference between neutron and gamma signals. 
Its low computational complexity makes the real-time dis- 
crimination of FPGA easy to achieve. The input and output 
signals of the trapezoidal shaper are processed separately in 
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Fig. 1. (a) Block schematic of the real-time digital pulse shape dis- 
criminator. (b) Real-time digital pulse shape discriminator. 


two parallel modules. One module is the n/y discrimination 
module, which identifies signal types by analyzing the char- 
acteristics of the shaped signals, and the other is the peak ex- 
traction module, which selects the maximum amplitude from 
the pulse data in the duration of the original signals. Conse- 
quently, with the identified signal types and amplitudes, the 
spectrum formation module forms gamma spectra, neutron 
response spectra, and PSD plots, which are updated in real- 
time after a new pulse signal arrives. Owing to the low power 


41 consumption of the discriminator and the small amount of 


spectrum data, only one USB port is required for power sup- 
ply and data transmission. 


B. Pulse-shaping technique 


In the radiation detection system, RC-CR circuits and a 
Sallen-Key filter were used to filter out noise from the ampli- 
fied nuclear signals to obtain an optimum signal-to-noise ratio 
(SNR). These filters remove the high-frequency components 
in the frequency domain, causing a change in the shape of the 
signals in the time domain. For example, RC-CR circuits and 
a Sallen-Key filter can process a signal into a Gaussian-like 
shape in both analog and digital ways [26-28]. In addition 
to filtering, the signals must be specifically shaped to accom- 
modate different application scenarios. The signal should be 
narrow at high count rates and have a flat top for amplitude 
analysis. In time analysis, steepening the leading edge of the 
signals makes the timing accurate. Using the digital pulse- 
shaping technique, a digital shaper can meet diverse require- 
ments and process signals into a variety of shapes, such as 
trapezoidal and triangular. 

A trapezoidal pulse-shaping algorithm was developed us- 
ing the z-transform method[29]. The z-transform method 
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comprises two main steps. The first step is to acquire the 
transfer function H(z) through the z-transformation of the in- 
put and output signals. The next step is to apply an inverse 
z-transform to obtain the equations for the input and output 
signals in the time domain. 

Common signal-processing algorithms, such as fast 
Fourier transform (FFT) and finite impulse response (FIR) 
filters, require numerous multiplication and addition opera- 
tions. Moreover, the internal DSP slices of the FPGA can 
conveniently perform basic multiplication and addition oper- 
ations. As shown in Fig. 2(a), the trapezoidal pulse-shaping 
algorithm derived from the z-transform method could be im- 
plemented as a shaper in FPGA. Because this algorithm re- 
quires high-precision multiplication and addition operations, 
the adders and multipliers should have large bit widths. This 
algorithm also requires additional operations to remove the 
baseline drift caused by noise accumulation [29]. 

Another method of achieving a digital trapezoidal shaper is 
the unfolding synthesis technique. Because of the limitations 
of the bit width, hardware resources, and real-time processing 
requirements, a trapezoidal shaper should be simple and effi- 
cient for FPGA. The unfolding technique was introduced to 
transform digital signals from an ADC into a unit impulse, 
which can approximately synthesize any pulse shape [17]. 
The nuclear signal amplitude increased rapidly to a maximum 
value and then decayed at a relatively slow rate. For Nal:T1 
crystals, the decay time of light output has both fast and slow 
components [30]. The current signal output from the PMT 
was converted into a voltage pulse in the resistive feedback 
preamplifier followed by an RC network. Fig. 3(a) shows the 
digitalized voltage pulse, which can be expressed as: 

Vin] =A (eF -e*), (1) 
where 7, and 72 represent the decay time constants of 
the digital signal in a bi-exponential shape, which includes 
the preamplifier time constant, scintillation fluorescence 
time constant, electronic characteristics of the anode, and 
preamplifier-anode coupling. A indicates the pulse ampli- 
tude and also controls the signal width. The interval between 
two adjacent sampling points n is determined by the sampling 
rate. Eq. (1) serves as a fitting model to obtain the decay time 
constants of the digital pulse. 

Because the output signal of the detector has a biexponen- 
tial shape, the first step is to deconvolute the biexponential 
signal into a unit impulse d[n]. The operational algorithm of 
a system that unfolds an exponential unit impulse response is 
given by the following equation[17]: 

öin] (2) 

There is also a double exponential impulse-shaping algo- 
rithm for obtaining a narrower 6 signal, which might result in 
better anti-pile-up applications[31]. All deconvolution meth- 
ods produce an approximate 6 signal that retains time and am- 
plitude information. The difference between these methods 
lies in the width of the 6 signal and complexity of the algo- 
rithm. Although Eq. (2) is derived for a single exponential 
pulse, good shaping results can be obtained by using Eq. (2) 
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twice. Fig. 3(c) and Fig. 3(d) show that the signal becomes 
significantly narrower and smaller after two deconvolutions. 

The next step was to synthesize a unit pulse using a trape- 
zoidal pulse. The operational algorithm for the synthesis sys- 
tem was deduced from the shape of the target waveform. A 
standard trapezoidal pulse can be described as 


(3) 


After differentiating the primary function y/n], the function 
value of each interval becomes constant, resulting in d/n] 


(4) 


The weighted and shifted unit impulses can be used to rep- 
resent the differentiation of y/n]. Therefore, an alternative 
expression for d/n] is as follows: 


an] = -(up] -un -8) -Ê 


ô[n] is the result of the difference operation of u[n], and 
the synthesis system based on recursive difference equations 
in the time domain[32, 33] is given as follows: 


A 


ul{n—J])- (5) 


‘(uln—m] 


pin] = = Oln] —6[n—k]) — 7 - ([n—m] —d[n—Jj]), ©) 
din] = pln] + dfn — 1], (7) 
yin] = dln] + yin — 1], (8) 


where k is the duration of the rising/falling edge of the 
trapezoidal pulse and |m — j| is the duration of the flat time 
of the trapezoidal pulse. |m — j| and k could be large enough 
so that the signal has enough time to rise to its maximum am- 
plitude, thus compensating for the ballistic deficit effect[34]. 
However, a long signal duration causes pile-up events when 
the count rate is high, which is unhelpful for n/y discrim- 
ination. In this study, we did not consider the pile-up situa- 
tion. Thus, we can only discard pile-up events by determining 
whether there are multiple peaks in the pulse duration. The 
tail pile-up pulses could be separated when only counting was 
considered. However, when a neutron pulse and gamma pulse 
are stacked together, it is difficult to separate them while re- 
taining their distinct original features. The challenge of sepa- 
rating the neutron and gamma signal pile-ups remains. There- 
fore, the parameters of the algorithm were set to prevent the 
shape of the trapezoidal signal width from exceeding the orig- 
inal signal width. Fig. 2(b) shows the block diagram of the 
trapezoidal shape algorithm. The parameter b is the ampli- 
tude factor that controls the digital gain. The parameters a, 
and ag are determined by the decay time constant 7 and sam- 
ple period AT: 
= e AT / Phy 
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Fig. 2. (a) Functional block diagram of the trapezoidal shaper using 
shaper using unfolding-synthesis technique. 


= e AT/72 i 


a2 (10) 


In Fig. 2, compared with the z-transform method, the trape- 
zoidal shaper derived from the idea of unfolding synthesis 
technique has fewer adders and multipliers. In reality, fewer 
operations of the trapezoidal shaper are attributed to the de- 
convolution of the digital pulse. Signal deconvolution shifts 
the target of the trapezoidal shape from a bi-exponential sig- 
nal to a 6 signal that retains the energy information of the 
original signal[35]. The implementation of signal deconvolu- 
tion reduces baseline shifts at high count rates and improves 
energy resolution[35, 36]. Hence, fewer operations do not 
necessarily result in poor performance. The recursive method 
for the trapezoidal shaper using Eqs. (2), (6), (7), and (8) is 
shown in Fig. 2(b). 


C. PSD algorithm based on a trapezoidal shaper 


The digitalized neutron and gamma signals sampled from 
the ADC are shown in Fig. 3(a). Neutron and gamma wave- 
forms were normalized to the same area for comparison with 
the same gamma-equivalent energy. We consider the rise time 
as the time when a signal crosses from a specific low-voltage 


the z-transform method. (b) Functional block diagram of the trapezoidal 
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threshold to a specific upper-voltage threshold (10% to 90%), 
and the fall time as 90% to 10% of the signal tail. The rise 
time of the gamma signal was 58 ns and its fall time was 470 
ns. The rise and fall times of the neutron signal were 46 ns 
and 413 ns. The two signals in Fig. 3(a) could be fitted using 
a bi-exponential function. The fast component of the gamma 
signal was 70 ns and its slow component was 200 ns. The 
fast and slow components of the neutron signal were 64 ns 
and 160 ns. It is evident that the neutron and gamma signals 
have different decay time constants. The difference between 
the two signals appears to be small when viewed by the hu- 
man eye. It can be presumed that the results obtained after 
integrating the signals directly will also be similar. The use 
of a trapezoidal shaper avoids the direct processing of similar 
original signals and analyzes the new features of the shaped 
signals. Because parameters a, and ag are set for gamma sig- 
nals, the trapezoidal shaper does not match the neutron sig- 
nals. Thus, the two types of signals formed different shapes, 
as shown in Fig. 3(b). The gamma signal is shaped into 
an approximate trapezoidal signal, whereas the output of the 
neutron signal has the characteristics of overshoot and under- 
shoot owing to the unmatched parameters a; and a2. Accord- 
ing to the output characteristics of the neutron signals from 
the trapezoidal shaper, the n/y and peak extraction modules 
cooperate in the following states: 
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(1) tı denotes the time at which the signal amplitude ex- 
ceeds the baseline value. During time tı, the maximum value 
of the original signals is determined through successive com- 
parisons and recorded. 

(2) The time tz from when the signal decays to the baseline 
until it returns to zero is recorded. 

(3) After completing the record of tz, the peak value is 
sent to the spectrum formation module to identify whether 
the peak value belongs to the neutron response spectrum or 
the gamma spectrum by calculating the correlation factor A: 


(11) 


A takes values between O and 1, and it characterizes the 
magnitude of the undershoot and the overshoot of the signal. 
When the neutron pulse signal width is relatively narrow, the 
signal decays quickly causing an undershoot and resulting in 
a long tə and short tı for peak extraction. Typically, À of the 
neutron signal is smaller than that of the gamma signal, which 
means that the neutron peak is on the left side of the gamma 
peak in the PSD plot. 

Low sampling rates affect the trapezoidal shaping results, 
but more importantly, they result in inaccurate measurements 
of tı and t2. Owing to the low-sampling-rate requirement of 
the shaper and the amplification of the difference between the 
neutron and gamma signals, this algorithm is characterized 
by its ability to work at low sampling rates compared to some 
mainstream algorithms[22]. 


N= tı /(tı + t2). 
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Fig. 3. (a) Averaged gamma-ray and neutron pulses from the NaIL 
detector. (b) Illustration of the PSD algorithm based on the trape- 
zoidal shaper. (c)First signal deconvolution of the input signal when 
aı is 0.9048. (d) Second signal deconvolution of the input signal 
when ay is 0.9048 and az is 0.7515. 


The neutron and gamma signals were processed in the same 
trapezoidal shaper and shaped into a nonstandard trapezoidal 
pulse. In Fig. 3(b), although the output of the gamma signal 
has a relatively flat top, the output of the neutron signal has an 
uneven top that could shift the neutron response peak to the 
right in energy analysis. Thus, using the parallel-processing 
capabilities of the FPGA, the peak extraction module matches 
the peaks of the original signals to the shaped signals. 
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II. EXPERIMENTAL SETUP 


A 50.8 mm x 50.8 mm ĉLi-enriched (95%) NaI:T1 scin- 
tillator with 1% Li doping was chosen because it has a high 
light output yield exceeding 30000/MeV[6]. A preliminary 
check of the detector response to gamma rays was performed 
using a !37Cs gamma source. The energy resolution was 
7.01% at 662 keV. The experiments were implemented us- 
ing a 30 mCi 74!Am-Be neutron source. Thermal neutrons 
were indirectly detected through °Li(n, a)?H reaction, which 
emitted 3H and a particles with a total energy of 4.79 MeV. 
Only part of the energy was deposited in the Nal:TI°Li crystal 
because of the luminous efficiency of a and 3H particles[37]. 
The cross section of the thermal neutron capture on Li was 
940 b. A major advantage of using °Li to capture neutrons 
is that the reaction has a high reaction energy, and a neu- 
tron response energy of 4.79 MeV can be easily distinguished 
from the gamma background. Additionally, discrimination 
between the neutron and gamma signals in the low-energy re- 
gion can be avoided using the nuclear reaction method. 

The count rates can differ when the distance between the 
241 Am-Be source and the detector varies. The count rate of 
gamma rays and neutrons was 1.8 kcps when the distance was 
set to 20 cm with a 6 mm thick lead shield, and 22 kcps with- 
out a lead shield. The 7*!Am-Be source was placed 20 cm 
away from the detector with a 6 mm thick lead shield for mod- 
erate reduction of the gamma background. Four arrangements 
were used to determine the appropriate PSD parameters and 
evaluate the pulse-shape discrimination ability of the discrim- 
inator. 

The trapezoidal shaper is completed in the first part. The 
waveforms of digitalized signals were extracted using an in- 
tegrated logic analyzer. By fitting the neutron and gamma 
signals with a bi-exponential function, we can obtain approx- 
imate decay time constants to build a trapezoidal shaper. Af- 
ter loading the original pulse data from the integrated logic 
analyzer into the MATLAB workspace, a shaper was built 
and simulated in Simulink of MATLAB. The decay time con- 
stants obtained by fitting with a biexponential function were 
float-point-type constants. To enhance the operation speed 
and save data memory in FPGA, all parameters were con- 
verted into a fixed-point type. The output results of each op- 
eration were rounded and saturated on an integer overflow be- 
cause of the limitation of fixed-point precision. This is why 
the precision of the algorithm in FPGA decreases as the oper- 
ations increase and the final shaping results are not as accurate 
as the theoretical calculation. An inaccurate operation output 
can cause a baseline shift and pulse shape distortion. The 
algorithm operations were most likely halted when an over- 
flow occurred. One way to solve these problems is to control 
the value and bit width of the operation results and parame- 
ters, such as a1, G2, and b. Because shaped signals can be 
observed in real time using an integrated logic analyzer, the 
parameters of the algorithm were manually adjusted after ob- 
servation. Each adder, multiplier, and accumulator in FPGA 
were also set to have different output bit-width limitations. 

The second part tested whether the algorithm could make 
signals distinctive and whether the spectrum-formation mod- 
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ule worked to form a PSD plot. The trapezoidal shaper 
was implemented in the discriminator to complete the 
neutron/gamma discrimination algorithm. The spectrum- 
formation module in FPGA was selected to record the PSD 
plot. The discriminator sampled the neutron and gamma sig- 
nals for some time to acquire PSD plots, which demonstrated 
different values for all counts. 

Neutron and gamma peaks are observed in the PSD plot. 
Assuming that the peaks follow a Gaussian distribution, we 
use a Figure of Merit (FoM), defined as 


S 
FWHM, + FWHM,’ 


FoM = (12) 
to determine the quality of discrimination. 

S denotes the distance between the neutron response and 
gamma peaks, while FW H M, and FW H M, represent the 
full widths at half maximum of the corresponding peaks. The 
FoM is high when there are few statistical fluctuations; the 
two peaks are significantly separated, and a high FoM indi- 
cates a good discrimination ability. According to the PSD 
plot, the neutron signals can be distinguished if \ is smaller 
than a threshold, Astandara. After Astandard Was determined 
and set in FPGA, the spectrum formation module was modi- 
fied to form two-pulse amplitude spectra. 

Subsequently, the determined parameters and modified 
modules in FPGA were used to acquire the final spectra in the 
third part. For instance, in the experiment at a 50 MHz sam- 
pling rate, a, = 0.9048, az = 0.7515, b = 1.5, k= 10, |m — j| 
= 20, and Astandard = 0.45. The signals were classified as 
neutron signals when A was < 0.45. In this section, two 4096 
channel spectra are uploaded from FPGA to a computer every 
two seconds. These parameters were set manually beforehand 
and could be fixed for the current experimental conditions. 
Changes in the source intensity, detector type, and detector lo- 
cation from the source might invalidate these parameters and 
deteriorate the discrimination quality. Therefore, the experi- 
mental steps described above must be repeated to obtain the 
best parameters when the detection conditions change. The 
fourth part is the shielding experiment conducted in the pre- 
vious two parts. A 1.5 mm thick cadmium plate was used to 
absorb neutrons, which allowed us to compare the discrimi- 
nation efficiency and ensure that only a few neutron signals 
were omitted. 


IV. RESULTS AND DISCUSSION 


As discussed previously, a PSD discriminator and its em- 
bedded algorithm were designed and implemented. With- 
out lead shielding, several input and output signals from the 
trapezoidal shaper were captured using an integrated logic an- 
alyzer. With 6 mm lead shielding, the probability of pile-up 
events occurring was reduced. The presence of massive low- 
energy gamma rays results in many low-amplitude pulses, as 
shown in Fig. 4(a). The signal shapes, with pulse widths 
shorter than 2 us, are shown in Fig. 4(b). A gamma sig- 
nal with a large amplitude was deliberately acquired. Af- 
ter trapezoidal shaping, the gamma signals feature a flat top, 
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whereas the neutron signals exhibit an overshoot at the top. 
The gamma and neutron signals exhibit small and severe un- 
dershoots, respectively. Thus, a signal with an amplitude of 
2836 can be identified as a neutron signal by the evident over- 
shoot and undershoot. In Fig. 4(a), the neutron signal has a 
shorter decay time and high amplitude and already has a mi- 
nor undershoot before pulse shaping. This can be eliminated 
by an analog circuit of pole-zero cancellation, which was not 
used in this experiment, considering the remaining original 
difference in the output signals of Nal:T] and °Li detectors. 
The fast decay of the neutron signal causes a baseline drift 
in the signal after the neutron signal. However, in Fig. 4(b), 
the neutron signal is restored to the baseline relatively quickly 
after pulse shaping. The signal after the neutron signal was 
less affected by the baseline drift. This demonstrates the abil- 
ity of the pulse-shaping technique to achieve pile-up decom- 
position. Fast shapers can shorten the rise and fall times of 
signals, reducing the occurrence of tail pile-up events. The 
shaped-signal width can be narrowed by adjusting the param- 
eters of the trapezoidal shaper or using a cusp shaper. For 
example, a flat top can be eliminated when |m — j| is set to 
zero, and the signal is processed into a triangular shape with 
a short signal width. 
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Fig. 4. (a) Amplified signals from the NaIL detector. (b) Pulse shap- 
ing of mixed neutron-gamma signals. 


The quality of discrimination using the pulse-shaping 
method with the Nal:TI°Li detector is shown in Fig. 5. These 
plots were uploaded directly to the FPGA. The distribution 
of counts with a correlation factor À approximates a Gaus- 
sian distribution, and there is a wide gap between the two 
peaks. The left peak is the neutron response peak, whereas 
the right peak is the gamma peak due to the large number 
of low-energy gamma rays. The left peak disappeared af- 
ter adding a 1.5 mm thick cadmium plate, proving that the 
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left peak was a neutron response peak. The right peak was 
taller in Fig. 5(b) than in Fig. 5(a), which also implies that 
cadmium has a high (n, y) reaction cross section[38] and ac- 
cordingly produces more gamma rays. Setting a threshold 
of 50 allowed an energy measurement from 57.195 keV to 
5.695 MeV, which covered the energy range of the neutron 
response. Additionally, setting a threshold taller than 0 helps 
filter out unconcerned low-amplitude noise. In one minute, 
the FoM is 2.88 at a 50 MHz sampling rate allowing for the 
determination of the sound parameters, which is sufficient for 
efficiently separating the neutron and gamma signals. This 
effective separation was attributed to the unique shape differ- 
ence between the gamma and neutron pulses amplified by the 
trapezoidal shaper. 


Moreover, the signals can be discriminated at lower sam- 
pling rates using suitable parameters. We set a, to 0.8825, 
az to 0.6997, and Astandara to 0.45 for a sampling rate of 40 
MHz. We set a; to 0.8187, a2 to 0.5647, and Astandara to 
0.4 for a sampling rate of 25 MHz. The width of the trape- 
zoidal pulse remained constant. An increase or decrease in 
the sampling rate did not affect the inherent decay time of the 
signal. From Eq. (9) and Eq. (10), we deduce that changes 
in the sampling rate affect a, and ag, and alter the results of 
the trapezoidal shaper. The fall time of the neutron signal was 
413 ns, and there were no more than 24 sampling points on 
the falling edge when the sampling rate was 50 MHz. Af- 
ter the sampling rate was reduced to 25 MHz, there were no 
more than 12 sampling points at the falling edge of the neu- 
tron signal. Similarly, there were fewer than three sampling 
points at the rising edge. Sometimes, only one sampling point 
was observed during the signal rise time. Thus, owing to the 
short duration of a nuclear signal, a low sampling rate affects 
the fidelity of digitalized signals. Although the trapezoidal 
shaper still operates normally at a 25 MHz sampling rate, the 
number of sampling points available for calculating A is re- 
duced by half at a 25 MHz sampling rate compared to a 50 
MHz sampling rate, which makes A less accurate. As shown 
in Fig. 5(e), the inaccuracy of À broadens the bottom of the 
peak and increases the distance between the two peaks. The 
dispersion of the counts makes the FWHM narrower. There- 
fore, the FoM values calculated using Eq. 12 are higher. The 
results show great discrimination ability at lower sampling 
rates when only counting is considered, but the quantization 
error caused at a 25 MHz sampling rate could deteriorate the 
energy measurements. Therefore, peak extraction should be 
applied to the original or output signals of the other shapers. 


The sampling rate is not the main factor affecting the FoM 
value but a constraint of the pulse discrimination method. 
In some methods, the scale of the sampling rate has certain 
requirements for computational accuracy. The adaptability 
of the method described above allows discrimination over a 
wide range of sampling rates. The parameter selection of the 
trapezoidal shaper largely determines the discrimination per- 
formance of the method. The decay time constants of the 
output signals of the different detectors were not the same. 
Therefore, this pulse-shaping method must fit the digitalized 
signals of the ADC in advance to obtain the approximate de- 
cay time constants. 
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Fig. 5. PSD plot obtained at a 50 MHz sampling rate (a) without 
and (b) with cadmium shielding. PSD plot obtained at a 40 MHz 
sampling rate (c) without and (d) with cadmium shielding. PSD 
plot obtained at a 25 MHz sampling rate (e) without and (f) with 
cadmium shielding. 


The neutrons emitted by the 744Am-Be source were not 
monoenergetic. The main reason why the neutron spectrum 
is continuous is the continuous energy of a particles produced 
by the a decay of 741 Am[39]. Additionally, the energy loss of 
q particles exists owing to the collisions of a particles with 
target nuclei and Coulomb repulsion interaction. Neverthe- 
less, the reaction products directly transitions from the ex- 
cited state to the ground state and the energy transferred to the 
reaction products is uniquely determined. The total energy of 
4.79 MeV emitted from °Li(n, œ)’H reaction is carried by 3H 
and a particles, when the kinetic energy of the neutrons is ne- 
glected. This energy is the response of thermal neutrons and 
thus should exhibit a single-peak distribution at the energy of 
4.79 MeV in the energy spectrum. For different scintillators, 
the luminous efficiencies of 3H and a particles were different 
for the different scintillators. In this experiment, the response 
energy deposited in the Nal:TI°Li detector was assumed to be 
less than 4.79 MeV. 

Before implementing the PSD algorithm, we plotted Fig. 
6(a) and Fig. 6(b) based on the mixed neutron-gamma ra- 
diation field for a measurement time of two minutes. The 
gamma-ray peak with the highest energy is the full en- 
ergy peak of 4.438 MeV emitted from the °Be(a, n)!?C* 
reaction[40]. Subsequently, a single-escape peak occurred at 
3.927 MeV and a double-escape peak occurred at 3.416 MeV. 
The ?4!Am-Be neutron source was surrounded by boron- 
containing polyethylene plates that were used to prevent neu- 
tron leakage. Therefore, neutrons are captured by hydrogen in 
boron-containing polyethylene, which leads to a gamma-ray 
peak at an energy of 2.23 MeV. Moreover, !°B is a suitable 
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Fig. 6. Response spectrum of mixed neutron and gamma signals not 
using the PSD algorithm (a) without and (b) with cadmium shield- 
ing. (c) Neutron response spectrum and (d) gamma spectrum ob- 
tained using the PSD algorithm without cadmium shielding. 


material for shielding neutrons from 1°B(n, a)" Li* reactions. 
The ‘Li in the excited state quickly deexcites to its ground 
state and emits gamma-rays with 478 keV energy. Relying 
on the high thermal neutron absorption cross section of Cd, 
the two spectra were combined to confirm the energy range 
of the neutron response. The area marked in red indicates the 
response energy range of the thermal neutrons that occur near 
the single-escape peak of !7C. 


After implementing the PSD algorithm, we plotted Fig. 
6(c) and Fig. 6(d), which show the separated spectra uploaded 
from the discriminator to the computer. The response peak for 
the reaction of thermal neutrons with °Li lies at an energy of 
3.9 MeV. The successful extraction of the neutron response 
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spectrum indicates that neutron signals with different ampli- 
tudes can be discriminated. Although the pulse width in- 
creased with the pulse amplitude, the pulse amplitude did not 
affect the intrinsic decay time constant. Accordingly, changes 
in the signal amplitude do not affect the recognition of the 
neutron signals. The separated gamma-energy spectra shown 
in Fig. 6(d) are similar to those in Fig. 6(b), but the counts 
will be higher in low-energy areas with cadmium shielding 
owing to the 1!3Cd(n, y) reaction. Then, we compare the red 
peak area in the four spectra, which is the energy range from 
3.75 MeV to 4.07 MeV. The total counts in the energy range 
with the gamma and neutron signals in Fig. 6(c) are 9987, 
and the counts reduce to 4903 when neutrons are shielded by 
cadmium. The peak areas of Fig. 6(c) and Fig. 6(d) are 5149 
and 4755, respectively. Their sum corresponds to the total 
counts of the spectrum in Fig. 6(a), which indicates that the 
pulse-shaping method does not cause significant count loss. 


V. CONCLUSION 


A portable n/y discriminator and its embedded algorithm 
were tested at low sampling rates. The results showed sat- 
isfactory PSD performances of the Nal:T1,°Li scintillator and 
discriminator. This method can also be used for other gamma- 
ray and neutron dual-detection scintillators, such as CLYC 
and CLLB. Unlike sophisticated computer-based PSD meth- 
ods, the low cost, low sampling rate, and real-time processing 
have significant advantages in practical applications. Future 
investigations into the application of this simple but robust 
discrimination method to the assembly of scintillator array 
and SiPM array are of particular interest. In radiation imag- 
ing applications, using n/y dual-detection scintillator arrays 
with this method can simultaneously and effectively obtain 
neutron and gamma images, thereby reducing the size and 
cost of the imaging instruments. 
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